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General-relativistic rotation laws in rotating fluid bodies 
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We formulate new general-relativistic extensions of Newtonian rotation laws for self-gravitating 
stationary fluids. They have been used to re-derive, in the first post-Newtonian approximation, 
the well known geometric dragging of frames. We derive two other general-relativistic weak-field 
effects within rotating tori: the recently discovered dynamic anti-dragging and a new effect that 
measures the deviation from the Keplerian motion and/or the contribution of the fluids selfgravity. 

One can use the rotation laws to study the uniqueness and the convergence of the post-Newtonian 
approximations, and the existence of the post-Newtonian limits. 
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I. INTRODUCTION 

Stationary Newtonian hydrodynamic configurations 
are characterized by a variety of rotation curves. The 
angular momentum per unit mass j can be any func¬ 
tion of r, where r is the distance from the rotation axis. 
Other restrictions arise from stability considerations 0. 
In contrast to that, for a long time the only known rota¬ 
tion law in general-relativistic hydrodynamics had been 
that with j as a linear function of the angular velocity. 
Recently Galeazzi, Yoshida and Eriguchi [2] have found a 
nonlinear angular velocity profile, that may approximate 
the Newtonian monomial rotation curves flo = w/r x 
in the nonrelativistic limit. In this paper we define 
general-relativistic rotation curves j = j(Q) that in the 
nonrelativistic limit exactly coincide with Qq = w/r x 
(0 < A < 2, A ^ 1). We are able to obtain the general- 
relativistic Keplerian rotation law that possesses the first 
post-Newtonian limit (1PN) and exactly encompasses the 
solution corresponding to the massless disk of dust in the 
Schwarzschild spacetime. 


II. HYDRODYNAMICAL EQUATIONS 

We recapitulate, following 0 , the equations of general- 
relativistic hydrodynamics. Einstein equations, with the 
signature of the metric (—, +, +, +), read 

R/m — 9n>^~2 = (1) 

where T is the stress-momentum tensor. The station¬ 
ary metric reads 


where p is the baryonic rest-mass density, h is the specific 
enthalpy, and p is the pressure. The 4-velocity u a is 
normalized, g a pu a u^ = —1. The coordinate (angular) 
velocity reads v = Qc> 0 , where f l = u^/u*. 

We assume a barotropic equation of state p = p(p). To 
be more concrete, one can take the polytropic equation of 
state p{p,S) = KT(S')p 7 , where S is the specific entropy 
of fluid. Then one has h(p,S) = K(S)^pjp 1 ~ 1 . The 
entropy is assumed to be constant. 

Define the square of the linear velocity 

V 2 =r 2 (n-^ye 2 ^~^ c \ (4) 

The potentials a , /3, v, and u> satisfy equations that have 
been found by Komatsu, Eriguchi and Hachisu 0. They 
constitute an overdetermined, but consistent, set of equa¬ 
tions. The general-relativistic Euler equations are solv¬ 
able, assuming an integrability condition - - that the an¬ 
gular momentum per unit mass, 


j = 


V 2 

(1 


VjU ’ 
c 2 J 


( 5 ) 


depends only on the angular velocity Q; j = j(Q). In 
such a case the Euler equations reduce to a general- 
relativistic integro-algcbraic Bernoulli equation, that em¬ 
bodies the hydrodynamic information carried by the con¬ 
tinuity equations = 0 and the baryonic mass con¬ 
servation = 0. It is given by the expression 
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2 J dfij(n) = c. 
( 6 ) 


ds 2 = —e~£ ( dx °) 2 + r 2 e~& (d<j) — ^ (r, z) dx 

+ e~£ ( dr 2 + dz 2 ) . (2) 

Here = ct is the rescaled time coordinate, and r, z, (j) 
are cylindrical coordinates. We assume axial symmetry 
and employ the stress-momentum tensor 

T a P = p{c 2 + h)u a uP + pg afi , (3) 


III. ROTATION LAWS 

The general-relativistic rotation law employed in the 
literature 0-0 has the form 

j(n) = A 2 (n c -n), (7) 

where A and Q c are parameters. In the Newtonian limit 
and large A one arrives at the rigid rotation, Q = Q c , 
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while for small A one gets the constant angular momen¬ 
tum per unit mass. A three-parameter expression for j 
is proposed in 

Below we define a new family of rotation laws, 


m = 


ud~ s n s 

1 _ J^yjl-SQl + S + * ’ 


( 8 ) 


where w, 8 , k and \P are parameters. The rotation curves 
(r, z) ought to be recovered from the equation 


w 1 - 5 ^ 6 V 2 

1-+ f " (fl-£)(l-£)’ ( } 

For 5^—1, the general-relativistic Bernoulli equation 
© acquires a simple algebraic form 


1+ 4W./ 1 _£ X 


i _ + * 

c c 


TTT3JS 


= c. 


( 10 ) 


We shall explain now the meaning and status of the four 
constants w, 8, k and \P. Assume that there exists the 
Newtonian limit (the zeroth order of the post-Newtonian 
expansion OPN) of the rotation law. This yields 


F! 0 = -^. (11) 

Thus w and <5 can be obtained from the Newtonian limit. 
Moreover, the constant w is any real number, while 8 
is nonpositive — due to the stability requirement QJ 
and satisfies the bounds — oo < 8 < 0 and 5^—1. These 
two constants can be given apriori within the given range 
of values. Let us remark at this point that the rotation 
law l|8|). and consequently the Newtonian rotation (EH), 
applies primarily to single rotating toroids and toroids 
rotating around black holes. In the case of rotating stars 
one would have to construct a special differentially rotat¬ 
ing law, with the aim to avoid singularity at the rotation 
axis. 

The two limiting cases 5 = 0 and 8 = — oo correspond 
to the constant angular momentum per unit mass (flo = 
w/r 2 ) and the rigid rotation (f2 = w), respectively. The 
Keplerian rotation is related to the choice of 8 = —1/3 
and w 2 = GM , where M is a mass [9j. The case with 
S = — 1 should be considered separately, but we expect 
that the reasoning will be similar. 

The values of n and \I/ are problematical. One possibil¬ 
ity to get them is to apply the PN expansion. The rota¬ 
tion law in the PN expansion scheme should not be given 
apriori, but is expected to build up — in the subsequent 
orders of c~ 2 — from the Newtonian rotation law. The 
Newtonian rotation curves are specified arbitrarily, but 
the next PN corrections should be defined uniquely. This 
is, however, a well known property of the post-Newtonian 


expansions, that they are non-unique. Damour, Jara- 
nowski and Schafer |8( demand that a test body rotating 
circularly in a Schwarzschild space-time satisfies exactly 
the Keplerian rotation law with fl 2 = GM/R 3 , where 
R is the areal radius. Inspired by this we impose a Fix¬ 
ing Condition (F-Condition thereafter) — that a rotat¬ 
ing infinitely thin disk made of dust in a Schwarzschild 
space-time satisfies exactly the Bernoulli equation and 
the Keplerian rotation law. 

Consider a rotating, infinitely thin and weightless disk 
of dust in the Schwarzschild geometry. This is a textbook 
knowledge that there exists a stationary solution — each 
particle of dust can move along a circular traj ectory o f 
a radius R with the angular velocity fl = yjGM / R?. 
We shall present this solution in conformal coordinates, 
using our formalism. The conformal Schwarzschild met¬ 
ric reads ds 2 = —<P 2 // 2 (dx °) 2 + / 4 (dr 2 + dz 2 + r 2 d(f > 2 ) , 
where $ = 1 — GM/(2c 2 Vr 2 + z 2 ) and / = 1 + 
GM/(2c 2 y/r 2 + z 2 ). The angular velocity is equal to 
the Keplerian velocity H 2 = GM/(\/r 2 + z 2 f e ) and 
R = y/r 2 + z 2 f 2 . The total energy per unit mass ’P van¬ 
ishes for a test dust. Let the disk lie on the z — 0 plane 
and assume the rotation law with 8 = —1/3 and k = 3: 

w 4/3 Q -1/3 y2 

l_^4/3Q2/3 “ft (!_£)■ (12) 

Here V 2 = fl 2 r 2 / 6 /$ 2 . Notice that h = 0; the enthalpy 
per unit mass vanishes for dust. This is a simple exercise 
to show that w = y/GM and SI 2 = GM/(r 3 / 6 ) solve both 
Eq. (fl2ll and the Bernoulli equation (flTH) : the constant in 
m equals to unity. 


IV. 1PN CORRECTIONS TO ANGULAR 
VELOCITY 

Taking into account the above, we shall prove 
that if k = (1 — 35)/(l + 8) + 0(c~ 2 ) and ’P = 
4co + (D(c~ 2 ), where co is the Newtonian hydro¬ 
dynamic energy per unit mass, then the exact 
solution satisfies the first post-Newtonian (1PN) 
equations. We shall use the formalism of Q and the 
rotation law © , and recover most of the results obtained 
in the 1PN approach employed in 0- Notice that if 
8 = —1/3, then k = 3 — one recovers the coefficient 
in front of zn 4 / 3 H 2 / 3 in <H2D that is required by the F- 
Condition. 

The 1PN approximation corresponds to the choice of 
metric exponents a = /3 = —v = —U with \U\ <C c 2 fic| . 
Define oj = r~ 2 A The spatial part of the metric 

ds 2 = — ^1 + + -j-'j (dx 0 ) 2 — 2c~ 3 A^dx 0 d(f> 

+ ^1 - (dr 2 + dz 2 + r 2 d4> 2 ) . (13) 

is conformally flat. 
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We split different quantities (p, p , h , U, and v l ) into 


their Newtonian (denoted by subscript ‘0’ and 1PN (de¬ 
noted by subscript T’) parts. E.g., for p , 4', and U 
this splitting reads 

P = Po + c~ 2 pi, 

(14a) 

fl = + c~ 2 v f, 

(14b) 

<F = 4«o + 0 (c- 2 ), 

(14c) 

U = U 0 + c~ 2 Ui. 

(14d) 

Notice that, up to the 1PN order, 


-d iP = dih 0 + c~ 2 dihi+0(c~ 4 ), 

(15) 


P 


where the 1PN correction hi to the specific enthalpy can 
be written as hi = pi. For the polytropic equation of 
state this gives hi = (7 — 1 ) hopi/po- 

Making use of the introduced above splitting of quanti¬ 
ties into Newtonian OPN and 1PN parts one can extract 
from Eq. d 6 j) the OPN- and lPN-level Bernoulli equations. 
The OPN equation reads 

ho + Uo- ^ + ^ n 2 r 2 = c 0 , (16) 

where Co is a constant that can be interpreted as the en¬ 
ergy per unit mass. At the Newtonian level this is sup¬ 
plemented by the Poisson equation for the gravitational 
potential 


A Uq = AirGpo: (17) 


where A denotes the flat Laplacian. The first correction 
vf to the angular velocity Q is obtained from the pertur¬ 
bation expansion of the rotation law © up to terms of 
the order c~ 2 . Assuming that \Po = 4 cq , one arrives at 


v 


<t> 

1 


2 2 , A^ 4fl 0 ^o 

l - ^ ° r r 2 {l-6) ~ T^T’ 


(18) 


where we applied Eqs. HID and m- 

Remember that in the Newtonian gauge imposed in 
the line element (fl3l) the geometric distance to the ro¬ 
tation axis is given by r = r( 1 — Uq/c 2 ) + 0(c~ 4 ). It 
is enlightening to write down the full expression for the 
angular velocity, up to the terms 0(c~ 4 ): 


fl — f In -4- — 


c 2 f2/(l-S) C 2( X _ , 5 ) 

Ad, 4 


Ho {Uo + Hpr 2 ) 


r 2 c 2 (1 — 6) c 2 ( 1 — 6) 

This expression reduces to 


floho- 


(19) 


H — S 2 o H—r — 


w 


f 2 /(i—<5) rc 2 (l — 6) 


flnh\ 


0 ' 1 0 > 


( 20 ) 


in the case of test fluids, at the symmetry plane 2 = 0 . 
For the dust, in the Schwarzschild geometry, we get 

= + ! = ( 21 ) 

the 1PN correction to fin is equal to Thus the 

FC condition is satisfied in the 1PN order. 

After these consideration we are able to interpret the 
meaning of various contributions to the 1PN angular ve¬ 
locity fl. The first term is simply the Newtonian rotation 
law rewritten as a function of the geometric distance, as 
given at the 1PN level of approximation, from the rota¬ 
tion axis. The second term in (fTTHl vanishes at the plane 
of symmetry, 2 = 0, for circular Keplerian motion of test 
fluids in the monopole potential —GM/R. Thus it is sen¬ 
sitive both to the contribution of the disk self-gravity at 
the plane 2 = 0 and the deviation from the strictly Kep¬ 
lerian motion. The third term is responsible for the geo¬ 
metric frame dragging. The last term represents the re¬ 
cently discovered dynamic anti-dragging effect; it agrees 
(for the monomial angular velocities Hq = r _ 2 ^ 1 _l 5 ^u;) 
with the result obtained earlier in [ll| . 

A comment on the term — ^y^jyHgr 2 , that has been 
missing in mi- The reason for this omission is following. 
There is a gauge freedom in choosing an integrability 
condition for the 1PN hydrodynamic equation; due to 
that the Bernoulli equation of the 1PN order is specified 
up to a function F(r). We assumed in El , in order to get 
the 1PN Bernoulli equation as in 0 , that F(r) = 0; but 
that is not consistent with the F-Condition. It appears 
that the right value is F(r) = — Vt^r 4 /{I + 5), which leads 
to the emergence of the term in question. 

The vectorial component A$ satisfies the following 
equation 


AA 0 - 2^^ = -16TrGr 2 p 0 n 0 . (22) 

The 1PN Bernoulli equation does not influence the 
1PN correction to the angular velocity. It has the form 

3 

ci = -hi -Ui- + 2r 2 (fl 0 ) 2 h 0 - h 2 0 

-Ah 0 U 0 - 2U 2 - 4 6 { ~+ 6) r 4 n 4 + F(r), (23) 

where Ci is a constant. In order to derive (1231) we again 
used 4>o = 4co- This result agrees with the 1PN calcula¬ 
tion of El U P to the term F(r). 

The 1PN potential correction U\ can be obtained from 

AZ7i = 4 t tG (pi + 2p 0 + p 0 (h 0 - 2U 0 + 2r 2 (n 0 ) 2 )) . 

_i 24) 

Equations 03) and EH) have been derived in |Tl| in 
the framework of 1PN approximation. They can be also 
obtained directly from the Einstein equations written for 
the metric ©, as derived e.g. in Q. Here we recall a 
version similar to that used in Q|; it turns out to be more 
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conveninent than the original form of Q. The relevant 
equations read 


Av = 47 t ^re 2a /' 


p ( c2 + h ^V 2 * /c 2 +2P 


+ 7 TT r 2 e 2 ^~ u)/c \io ■ Vw - + u)-Vu 

2 c 4 c 1 


and 


A 


+ -d r ) w = -16 t r^e 2 “ /c2 p(c 2 + fe) ° 

r ) c z 1 — V z /c 


+ -V(t/- 3/3) • Vw, 


where V denotes the “flat” gradient operator. The re¬ 
maining Einstein equations yield corrections of higher 
orders. 

In summary, we have shown that — for —oo < 
w < oo and —oo < S < 0, 8 ^ — 1 — the choice 
n = (3 — <5)/(l + <5) + 0(c~ 2 ) and = 4co + 0(c~ 2 ) 

in the rotation law <[9|) guarantees that if there 
exists an exact solution analytic in powers of c -2 , 
then it satisfies the OPN and 1PN approximating 
equations. 

One easily finds out that the rotation law (0) satisfies 
the generalized Rayleigh criterion 0 ] for stability ^ < 
0 up to 1PN order, assuming that S is strictly negative. 

Comments on the 1PN corrections to the angular ve¬ 
locity. In the following considerations we assume w > 0, 
which means fly > 0 , but the reasoning is symmetric un¬ 
der the parity operation w —> —w. The specific enthalpy 
h > 0 is nonnegative, thence — is nonpositive 

the discovered in 0 instantaneous 1PN dynamic reac¬ 
tion slows the motion: it “anti-draggs” a system. In 
contrast to that, the well known geometric term with A^ 
is positive p~T| , and the contribution to the angu¬ 

lar velocity is positive it pushes a rotating fluid body 
forward. Thus the two terms in (1181) counteract. 

Dust is special - the specific enthalpy ho vanishes, 
hence dust test bodies are exposed only to the geometric 
effect — the frame dragging. Even more special is the 
rigid (uniform) rotation — the correction terms vf are 
proportional to 1/(1 — 5) and they vanish, because now 
S = —oo. Uniformly rotating disks are already known to 
minimize the total mass-energy for a given baryon num¬ 
ber and total angular momentum [l3|- The vanishing 
of the 1PN correction vf is their another distinguishing 
feature. 

It follows from our discussion that assuming the F- 
Condition, one has three free parameters: w, 5 and 4'; 


the parameter n is a given function of S. The full system 
of Einstein-Bernoulli equations can be solved numerically 
within this class of data and the resulting solutions are 
expected to possess OPN and 1PN limits. 

V. CONCLUDING REMARKS 


We write down the general-relativistic rotation laws, 
recover the well known geometric dragging of frames 
and derive a full form of the two other weak-fields ef¬ 
fects, including the dynamic anti-dragging effect of 0 - 
The latter can be robust according to the numerics of 


1 lj, but the ultimate conclusion requires a fully general- 


relativistic treatment, that is the use of the new rota¬ 
tion laws. The frame dragging occurs — thr oug h the 
Bardeen-Petterson effect 0 ] - in some AGN’s [l7j |. The 
two other effects can lead to its observable modifications 
in black hole systems with heavy disks. 

In the weak field approximation of general relativity 
the angular velocity of toroids depends primarily on the 
distance from the rotation axis — as in the Newtonian 
hydrodynamics - but the weak fields contributions make 
the rotation curve dependent on the height above the 
symmetry plane of a toroid. 

The new rotation laws would allow the investigation 
of self-gravitating fluid bodies in the regime of strong 
gravity for general-relativistic versions of Newtonian ro¬ 
tation curves. In particular, they can be used in order 
to describe stationary heavy disks in tight accretion sys¬ 
tems with central black holes. These highly relativistic 
systems can be created in the merger of compact bina¬ 
ries consisting of pairs of black holes and neutron stars 
00 , but they might exist in some active galactic nu¬ 
clei. 

The new general-relativistic rotation laws can be ap¬ 
plied to the study of various open problems in the post- 
Newtonian perturbation scheme of general-relativistic 
hydrodynamics. We demonstrate in this paper that an 
adaptation of the condition used in Q ensures uniqueness 
up the 1PN order. Further applications include the in¬ 
vestigation of convergence of the post-Newtonian pertu- 
bation scheme, as well as the existence of the Newtonian 
and post-Newtonian limits of solutions. 
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